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Local electric current correlation function in 
an exponentially decaying magnetic field 

N. SadooghEl and F. TaghinavaJll 
Department of Physics, Sharif University of Technology, P.O. Box 11155-9161, Tehran-Iran 

The effect of an exponentially decaying magnetic field on the dynamics of Dirac fermions in 3 + 1 
dimensions is explored. The spatially decaying magnetic field is assumed to be aligned in the third 
direction, and is defined by B(x) = B{x)ez, with B{x) — Bqc'^^'^^ . Here, ,f is a dimensionless 
damping factor and (.b = (eBo)^^''^ is the magnetic length. As it turns out, the energy spectrum 
of fermions in this inhomogeneous magnetic field can be analytically determined using the Ritus 
method. Assuming the magnetic field to be strong, the chiral condensate and the local electric 
current correlation function are computed in the lowest Landau level (LLL) approximation and the 
results are compared with those arising from a strong homogeneous magnetic field. Although the 
constant magnetic field -Bo can be reproduced by taking the limit of ^ — >■ and/or x — ^ from B{x), 
these limits turn out to be singular, especially once the quantum corrections are taken into account. 

PACS numbers: 11.15.-q, 11.30.Qc, 11.30.Rd, 12.20.-m, 12.20.Ds 



I. INTRODUCTION 

The effects of magnetic fields on systems containing 
relativistic fermions are subject of intense theoretical 
studies over the past two decades. These effects in- 
clude a wide range, from condensed matter physics [1[ 
to high energy physics [2| and cosmology [^. Theo- 
retical studies deal in the most of these works with 
the idealized limit of constant and homogeneous mag- 
netic fields. However, this limit is only reliable as 
long as the scale of field variation is much larger than 
the Compton wavelength of the fermionic system. As 
it is shown in [^, strong and homogeneous magnetic 
fields lead to the formation of fermion bound states 
even in the weakest attractive interaction between 
the fermions. This phenomenon, known as magnetic 
catalysis of dynamical chiral symmetry breaking [5| , is 
essentially based on a dimensional reduction from D 
to D — 2 dimensions in the presence of strong homoge- 
neous magnetic fields in the regime of LLL dominance. 
However, as it is demonstrated recently in [3] , massless 
Dirac electrons can also be confined by inhomogeneous 
magnetic fields in graphene, which is a single layer of 
carbon atoms in a honeycomb lattice. Electrons in 
graphene are described by a massless two-dimensional, 
time-independent relativistic Dirac equation. As it is 
shown in [6|, magnetic quantum dots can be formed, 
if one neglects the effects of electron spin and solves 
the Dirac equation in the presence of a square well 
magnetic barrier of width 2d, B(x) = B{x)ez, di- 
rected perpendicular to the graphene x — y plane, with 
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B{x) — Bo9{(P — x^), and 9 the Heaviside step func- 
tion. A circularly symmetric magnetic quantum dot 
arises then by a radially inhomogeneous magnetic field 
B = B{r)e^ (see 0,0] for more details). In ^, the 
bound state solutions and the spectra of graphene ex- 
citations in the presence of an inverse radial magnetic 
field, with B{r) ^ -, is studied. Analytical solutions 
for the Dirac equation in the presence of double or 
multiple magnetic barriers are presented in [9|. The 
bound spectra of electrons in a magnetic barrier with 
hyperbolic profile is studied in [lOJ. 

In the present paper, we will focus, in particular, 
on the solution of Dirac equation in the presence of 
an exponentially decaying magnetic field. Analytical 
solution of a quasi-two-dimensional Dirac equation in- 
cluding the spin of a single electron in a magnetic field 
of constant direction with arbitrarily strong exponen- 
tially depending variation perpendicular to the field 
direction is firstly presented in llj . In [12| , the prob- 
ability density and current distributions of Dirac elec- 
trons in graphene in the presence of an exponentially 
decaying magnetic field are determined. The forma- 
tion of chiral condensate induced by exponentially de- 
caying magnetic fields, in 2 -I- 1 dimensions, is studied 
in |l3| . The vacuum polarization tensor of 3-1-1 dimen- 
sional scalar Quantum Electrodynamics in the pres- 
ence of an inhomogeneous background magnetic field, 
arising from a general plane wave basis for the under- 
lying gauge field, is computed recently in [IJ]. In order 
to demonstrate the nonlocal features of quantum field 
theory in the presence of inhomogeneous background 
magnetic fields, the vacuum polarization tensor is de- 
fined to be local. This turns out to be a useful task 
providing local information about the nonlocal nature 
of fluctuation- induced processes [ij]. In the present 
paper, following the same idea as in [lj|, we will de- 



fine a local electric current correlation function (local 
electric susceptibility), and compute it in the presence 
of a strong and exponentially decaying magnetic field 
in the LLL approximation. Eventually, we will com- 
pare the result with the electric current correlation 
function arising from a strong homogeneous magnetic 
field and discuss the effects of inhomogeneity of the 
background magnetic field. 

The organization of the paper is as follows: In Sec. 
Ill Al we will briefly review the Ritus eigenfunction 
method [l5| by solving the Dirac equation in the pres- 
ence of a constant and homogeneous magnetic field. In 
Sec. IIIBl we will introduce an exponentially decaying 
magnetic field B(a;) = B{x)ez, which is assumed to 
be aligned in the third direction, with B{x) given by 
B{x) — _Boe~^^/^^. Here, ^ is a dimensionless damp- 
ing parameter and £b ~ (e_Bo)^^" is the magnetic 
length. Using the Ritus method, the energy eigenval- 
ues and eigenfunctions of a 3 -I- 1 dimensional Dirac 
equation in the presence of B(x) will be determined. 
In Sec. mil we will determine the chiral condensate 
{ip'>p)j using the fermion propagator in the presence 
of strong homogeneous (Sec. IIII Ap and exponentially 
decaying magnetic fields (Sec. IIIISp in the LLL ap- 
proximation. We will show that whereas (ipip) in a 
strong homogeneous magnetic field is constant, (ipip) 
in an exponentially decaying magnetic field depends 
nontrivially on a dimensionless variable m, defined by 
u = ^e~^^/^^. In Sec. IIVI the local electric cur- 
rent correlation function (local electric susceptibility) 
^(') will be computed in the presence of strong ho- 
mogeneous (Sec. IIV A[) and exponentially decaying 
magnetic fields (Sec. IIVB[) in the LLL approxima- 
tion. Here, i = 1,2,3 denotes the three spatial direc- 
tions. As it turns out, in the regime of LLL dominance 
the electric susceptibility in the perpendicular direc- 
tion relative to the external magnetic field (i.e. for 
i = 1,2) vanishes in both homogeneous and inhomo- 
geneous magnetic fields. The longitudinal component 
of x^*^ denoted by x", turns out to be constant in a 
homogeneous magnetic field and depend nontrivially 
on the dimensionless variable u in an exponentially de- 
caying magnetic field. In Sec. [Vj we will compare the 
results for the chiral condensate and the local electric 
current correlation function arising from homogeneous 
and inhomogeneous magnetic fields and show that the 
limits ^ — >■ as well as x — ^ are singular. In other 
words, although in the limits ^ — ^ and/or a; — )■ 0, we 
get B{x) — )■ Bq, but the values of (f/'V') and x" arising 
from an exponentially decaying magnetic field in these 
limits are not the same as {'iptjj) and x arising from a 
homogeneous magnetic field. This demonstrates the 
nonlocal features of the quantum vacuum in a strong 
and spatially decaying magnetic field. As is indicated 



in Sec. IVI[ which is devoted to our concluding re- 
marks, the results of this paper are relevant for high 
energy physics as well as condensed matter physics. 



II. ELECTRONS IN EXTERNAL MAGNETIC 
FIELDS 

Ritus eigenfunction method Il5| is a powerful tool to 
study the dynamics of electrons in external electro- 
magnetic fields. It is recently used in [16] to solve the 
Dirac equation in the presence of a uniform magnetic 
field in 2+1 dimensions, a uniform electric field in 1-fl 
dimensions and an exponentially decaying magnetic 



field in 2 -t- 1 dimensions (see also 13|). In IT], |l8|, 
the Ritus method is used to solve the 3-1-1 dimensional 
Dirac equation for massless fermions in the presence 
of a constant magnetic field. In Sec. Ill Al we will 
briefly review the results from 18] and will present 
the solution of the Dirac equation for a single massive 
electron in the presence of a constant magnetic field in 
3 + 1 dimensions. This will fix our notations. In Sec. 
IIIBl we will use the same method and will solve the 
Dirac equation for a massive electron in an exponen- 
tially decaying magnetic field. Our results coincide 



with the result presented in [13|, although our nota- 
tions are slightly different. We will also present the 
final form of the fermion propagator in the presence 
of a uniform and an exponentially decaying magnetic 
field in Sees. Ill Al and IIIBl respectively. The results 
presented in this section will then be used in Sec. IIIII 
to determine the chiral condensate, and in Sec. IIVI 
the local electric current correlation function in the 
presence of uniform and non-uniform magnetic fields. 



A. Electrons in constant magnetic fields 

Let us start with the Dirac equation of a single elec- 
tron in 3 + 1 dimensional Minkowski space in a back- 
ground electromagnetic field 



(7-n-m)V' = 0, 



(ILl) 



where H^ = i9^ — eA^ and A^ is the electromag- 
netic potential. Here, m and e are the electron's mass 
and electric charge, respectively. The aim is to solve 
pi.ip and to eventually determine the fermion prop- 
agator in the presence of background constant mag- 
netic field B = BqBz, aligned in the third direction. 
Such a constant magnetic field, is built, for instance, 
by choosing the gauge field A^ in the Landau gauge 
as Ag = (0,0, Boa;, 0). To solve (jILip . we use, as in 
16|,ll7|, the Ansatz ip — EpUp. Here, Ep is a diagonal 



matrix satisfying 



(7.n)Ep = Ep (7.p), 



(II.2) 



and Up is a free spinor, that describes an electron with 
momentum^ and satisfies {•j-p—m)up — 0. The Ritus 
eigenfunction Ep and the Ritus momentum p in (jll.2|) 
are unknown and shall be determined in what follows. 
The matrix Ep is determined by solving the eigenvalue 
equation 



(7 • n)^ Ep = p'Ep, 



(II.3) 



that arises directly from (jll.2|) . Using (7-11)^ = E^ + 
^cr^^Ff,^, with a'"' = |[7'",7''] and the field strength 
tensor F^i, = d^iA^ — di^Af^, with non- vanishing ele- 
ments F12 = — F21 = Bq, and choosing, without loss 
of generality, the Dirac 7-matrices as 7° = 0-3 ® 0-3, 
and 7* = ias (8) cr^, where ai,i = 1,2,3 are the Pauli 
matrices, (jll.3[) reads 



(n^ + cr^^eSo) Ep = fl 



(II.4) 



Here, a^'^ = 47^7^ = I (g) 0-3 with I a 2 x 2 identity 
matrix. To solve (jll.4[) . we use the fact that in the 
Landau gauge, the operator (7 • 11) commutes with 
Ho, Uj,j = 2, 3 and -P = - (7 • Uf + H^ - Uj, and 
therefore has simultaneous eigenfunctions with these 
operators. The eigenvalues of these operators, defined 
by the eigenvalue equations 



HoEp 
VEp = pEp, 



PoEp, IIjEp ^pjEp, j = 2,3, and 



(11.5) 



label the solutions of Dirac equation (III. II) in the back- 
ground magnetic field. Using (jll.3|) . the definition of 
V and (jII.Sp , it can be shown that the four-momentum 
p satisfies p^ = Pq—p^ —p, and T^Ep = p Ep is therefore 
given by 

{dl - {p2 - eBoxf + CT^^eBo) Ep = -p Ep. (II.6) 



Before solving (III. 61) , let us determine the components 
of Ritus momentum p, which turns out to play an 
important role in this method. To do this, we consider 
pi.2p . and use the Ansatz 



E„ 



/ , ^<jEp^„, 



(11.7) 



CT=± 



with the projectors 17 1 

A,^ = diag ((5^,1, (5^,-1, (5^,1, 5o-,-i) • (H.i 



Plugging first Ep from pi.7p in the left hand side 
(l.h.s) of (|II.2p . and then comparing the resulting ex- 
pression with the expression arising from the right 
hand side (r.h.s.) of (|II.2p . where again Ep from pi.7p 



and the Ansatz p — (po,0,p2,P3) is used,"'^ we ar- 
rive at (|II.6p . with p given by p = p^- This will fix 

Let us now turn back to the solution of (|IL6p . 
whose form suggests a plane wave Ansatz in the 0, 2 
and 3 directions. 



i?P,±i(x) = e-*^-/p-(x) 



(11.9) 



wher e x = (t, x,y,z) and p = (po, 0,p2,P3). Plugging 
(HlJt in pTf)) . and using A± = i (l ± «7^7^) for the 
projectors defined in (jILSp . the general solution for Ep 
is given by 



Ep(x) = e-'^Po*-P^y-P'''^Pp{x), 



(ILIO) 



where 



Pp{x) 



2{[fpi^) + fpi^)]+nV[f^{^) ~ /p-(x)]}.(ILll) 

To determine f^{x), it is enough to replace Ep(a;) 
from pLlOp in pL6p to arrive at 

(a^ - (p2 - eBoxf ± eBo) f^{x) = -p/p±(x).(IL12) 

Renaming the discrete quantum numbers y/p in the 
four-momentum p by .^/p = sgn(ei?o)-\/2|ei3o|p, and 
choosing a new coordinate C = ^ (2; — ^5^2)1 with Hb 
the magnetic length defined by Hb = \gBq\~'^/'^, the 
differential equation in pL12p can be reformulated in 
the form of the differential equation of harmonic oscil- 
lator, whose center lies in xq = P2/eBo and oscillates 
with the cyclotron frequency ujc = IbBq [16| . 
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As it turns out, the solution to the above equation 
can be given in terms of parabolic cylinder function 



0„„(a;) = a„„exp ( -^ j Hn^ f — 



with 



1 



where n„ is, in general, given by 

(T 1 

na=p+ -sgn(eBo) - -. 



(11.14) 



(IL15) 



^ The form p^ = Pq — P3 — p suggests that pi = 0. Moreover, 
we have fixed p; = pi,i = 0,3, where according to llII.SI l 
Pi, i = 0, 3 are the eigenvalue of Ili, i = 0, 3, respectively. 

^ Note that the general solution for p2 is either +y/p or —y/p. 
Here, we have chosen p2 = — \/p to have the same notation 
as in 1181 . 



Here, p labels the energy (Landau) levels. In the rest 
of this paper, we will work with ei?o > 0. According 
to pi.lSI) . for sgn(ei3o) > 0, the condition that n„ >0 
for p = 0, fixes the spin orientation of the electrons in 
the LLL to be positive (cr = +1) 
be given by 71+ = p for all values of p 



In this case n+ will 



17|. For this 



specific choice, the final result for f^{x) is the same 
as the result presented in 18], 



/p"(^) 



(f>p (x 



yp-i 



4p2) , 

^ - 4p2) 



P = 0,l,2,---, 
p=l,2,3,---. 
(11.16) 



Hence, for cBq > 0, only the positive spin solution 
/g (x) contributes in the LLL, characterized by p = 0. 
In other words, fgix) = fli{x) is undefined and is to 
be neglected, whenever it appears in a computation. 
The energy dispersion relation for an electron in 
the presence of a constant magnetic field is determined 
by the Ritus momentum 



Pp 



and is given by 



= fpo,0,-v/2^ 






2eBoP 



Pi 



0P,P3 



P-0,1,2, 



(11.17) 



(H.18) 



The Ritus eigenfunction Ep from (jII.10[) can be used 
to determine the fermion propagator in the presence 
of a uniform magnetic field. To do this, the Ansatz 
■ip = EpMp for the electron is to be generalized to a sys- 
tem including particles and antiparticles. Defining the 
associated creation and annihilation operators, and 
following the standard steps to determine the propa- 
gator 19[, the fermion propagator in the presence of 
a constant magnetic field is given by (see also [18,]) 



G{x,x') = WxMx')) 

Y^ /ppe'^P- (*-*') Pp(x) 



p=0' 



{l-Pp- m) 



Pp{x'), 
(H.19) 



where Vp = ^^j^P-, Pp{x) is given in pLTTj) with 
fp{x) from pi.lGp . and the Ritus momentum pp is 
given in pi.l8|) . The same expression appears also 
in [18| . Since we are working with eBo > 0, and 
the spin orientation of the electrons are already fixed 
in the LLL to be positive, no spin projector appears 
in pi.l9[) . as in the fermion propagators presented in 
n\ . To understand this, let us notice again that the 
facts that LLL includes only the positive spin solution, 
and that negative spin solutions contribute only to 
the higher Landau levels, are explicitly implemented 
in the choice p = 0, 1, 2, • • • for positive spin solution 



fp{x), and p = 1, 2, 3, • • • for negative spin solution 
f-{x) = /+_i(a;) in dlLiel) . In Sec. HHl we wih use 
the above method to determine the energy dispersion 
relation and the propagator for electrons in the pres- 
ence of exponentially decaying magnetic fields. 



B. Electrons in exponentially decaying magnetic 
fields 

In this section, the Ritus eigenfunction method will be 
used to determine the energy levels of a single fermion 
in the presence of an exponentially decaying magnetic 
field B(a;) = i?oe~"^e^, aligned in the third direc- 
tion. Here, Bq is the magnetic field at a: = 0, a is a 
dimensionful damping parameter. Later, we will set 
a = ^^/cBq, where ^ is a given dimensionless damp- 
ing parameter and eBo > 0- The dynamics of this 
electron in such an inhomogeneous magnetic field is 
described by the Dirac equation (jll.l[) . To solve it, 
we fix the gauge, in contrast to the previous case, by 
Ap{x) = (0,0,-^(e-"--l),0), as in [13. With 
this choice of Af^ , the case of constant magnetic field 
will be recovered by taking the limit a — > (or equiv- 
alently ^ ^- 0) in the classical level. Using the Ansatz 
ip = EpUp, with Ep satisfying (III. 21) . and following the 
same steps leading from (jll.2|) to (jll.6p . we arrive at 



dl 



idy 



eBo 



(e-""-l) +aeBoe- 



-pEp, 



E„ 



(H.20) 



where a — ±1 are the eigenvalues of third Pauli ma- 
trix, (73 — diag(l,— 1). Using in analogy to (|II.7p . a 
plane wave Ansatz for E„,'^ 



Ep{p) = J2 ^-EpAp)^ 



{11.21) 



<7=± 



with A„ defined in dlLSl) . and Ep^^{p) = e'P-^F^^{u), 
we arrive first at 



52 



d u^ 



du^ du 4 



P2 , cr \ Ap-Pl) 

U H ;; 



Fl{u)^Q. (11.22) 



In the above equations, p = (t,u,y, z) where w is a 
dimensionless variable defined by u = ^ei3oe~"^. 
Moreover, P2 ~ P2 + <^Bo/a, and {x,p) are given as 



Note that the plane wave Ansatz is justified by ['P,id'}] 



[V, ido] = 0, where V is defined by V ■ 
as in Sec. Ill Al 



'{7-n)^ 



-n^ 



n^ 



in the previous section [see below pL9|) ]. Comparing 
(jll.22[) with the differential equation 






d_ 
du 

(2n- 



u 

T 



satisfied by 



{n + ky. 



$,^;(u)=0, (11.23) 



-"/2efe/2£^(u), (11.24) 



where C^{u) is the associated Laguerre polynomial, 
the quantum numbers k and Tig- are given by 



k ^ 



pI^p, 



and 



1 , P2 - Vp2 -P 



(11.25) 



Note that the indices k and rio- in £jij are to be posi- 
tive integers. This implies the following quantization 
for p2 and p 



L|J 



s, 



and 



_ 2/ 2 2\ 

p = a [s — r ), 



leading to fc = 2r, with r > 0, and 



(T 1 

"^=2^2 



(11.26) 



(11.27) 



In (jll.26|) . [a\ is the greatest integer less than or equal 



11 



to a. As it is shown in 

ated with certain constant length £o 

the relation 



13j, s = L-P2/aJ is associ- 
through 



/7b^- 



1^1 = 1^1 



(11.28) 



where xq is a fixed length. At this stage, let us com- 
pare, Tier from (jll.27|) with n^- appearing in (JII.ISP for 
cBq > 0.^ As it turns out, in the case of exponentially 
decaying magnetic field, s — r labels the energy levels. 
In analogy to the case of constant magnetic field, to 
keep rirj from (|II.27I) positive, the electron spin orien- 
tation a in the lowest energy level, characterized by 
s ~ r ~ 0, is fixed to be positive (cr = -1-1). In this 
case, we have n = n^ = s ~ r. What concerns the 
quantum numbers r and s, for r > 0, we get s > r, 
that guarantees n — s — r > 0. In analogy to the 
solutions pi.lSp . for the constant magnetic fields, we 
arrive therefore at 

F^^u) EE Fll (u) = N,,^,e-^u^Cl^iu), (11.29) 



Comparing to llII.lSl l. ria from ||II.27|I . is in general given by 



with n = s — r = 0,l,2,--- , for positive spin solutions 
and 

F?:i,{u) EE F?:Ju) = iV„_i^,e-ty'-£2;_^(„)^ (n.30) 

with ri = s— r = l,2,3,---, for negative spin solutions. 
Using the orthonormality relations of associated Leg- 
endre polynomials C^ , the normalization factor Nn^r 
is given by 



Nn,r = 



{n + 2ry.' 



Using finally the projection matrices A± = ^(1 ^ 
17^7^), the solution for E„ can be brought in the form 
similar to the Ritus eigenfunction (III.10[) and pi.ll[) 
for constant magnetic field. 



Ep{p) = e-'(?'o*-?'^a-P^^)p,f (w). 



(11.31) 



with 



l{[F^:iu) + F^L,iu)]+^l'l'[F?:iu)~F?:_,iu)]}. 

(11.32) 

Using the orthonormality conditions for the associ- 
ated Laguerre polynomials, C^{u), presented in App. 
El it is straightforward to derive the closure and the 
orthonormality relations for Ep, 



oo oo 



r — 1 s—r 

as well as 



d2p|| 

(2^ 



E.p{p)Epip') = 6ip-p'), (11.33) 



VpEp{p)Ep'{p)^a-^S\p\\ -p\,)6s,s'K,,. (11.34) 



In the above relations, a = ^y/eBo, ^p{p) = 
7oE^(/o)7o, P|| = (poiPa) and S{p-p') = 5{t~t')6{u- 
u')5{y — y')5{z — z'), 'Dp = dtdudydz, and 



n^,^, =A+SryS,^r + {l-6s,r) 



0(r' - r)(-l)'-'-^/lS±ll£iiEl[ 




X A 



,(11.35) 



jSgn(eBo) 



with A± given in pi.Sp . These relations are similar 
to the results presented in [13|. Note that the term 
proportional to 5r,s isolates the contribution of the 
lowest energy level, characterized by s = r. This term 
is also proportional to A_|_, implying that the spin 



orientation in the lowest energy level for cBq > is 
positive. What concerns the contribution of higher 
energy levels, we had to distinguish between r = r' , 
r > r' and r < r' cases. This is because plugging Ep 
from (|IL3T|) and (|lL32l) in the l.h.s. of pOi)) . leads 
in general to the orthonormality relations similar to 
(jA.ip in App. [3 where r and r' are not necessarily 
equal. In App. \^ the orthonormality relation (jA.l[) 
is determined separately for r = r' , r > r' and r < r' . 
Using these results, it is straightforward to derive the 
remaining terms in (|II.35p , arising from higher energy 
levels. 

To determine the Ritus momentum p, we insert 
the above solution pOTI) and (llL32|) for Ep in (|IL2|) . 
After some straightforward algebraic manipulations, 
where we use the recursion relations of the Laguerre 
polynomials C'^{u), we arrive at 



(7 • n)Ep = Ep[7V + -/^ay/n{n + 2r) - j^ps]. (11.36) 

Here, n — s — r and a = ^^/cBq. Compar- 
ing the r.h.s. of (III. 361) with pi.2p . we get p — 
Ipo, 0, —OL^Jn (n + 2r),pz ) , which can alternatively 
be given in terms of (s, r) and ^ as 



Pk = 



{pa,0,-^\/eBoK,p3 



(11.37) 



with K 



^ — r^. Using (III.37p . the energy disper- 
sion relation of a single electron in an exponentially 
decaying magnetic field is given by 



£^ = JeeBoK 



pI + w?, 



(11.38) 



which is to be compared with the energy dispersion 
relation of an electron in the presence of a constant 
magnetic field from (jILlSp . Let us notice, that al- 
though n — s ~ r labels the energy levels, in contrast 
to the case of constant magnetic field, k — s'^ — r'^ 
determines the ordering of energy levels. Because of 
the above mentioned conditions for the integer quan- 
tum numbers r > 0, s > r (or 71 > 0), some value of 
K are not allowed. The simplest way to determine k, 
is to write it in terms of n and r as k = n(n + 2r). 
Choosing r — 1, 2, • • • , then k turns out to be k = 
for 71 = (or all values of r = s), k = 3, 5, 7, • • • , l + 2r 
for n = 1, K = 8, 12, 16, • • • , 2(2 -I- 2r) for n = 2, etc. 
In Sees. Illll and llVl the summation over energy levels 
in an exponentially decaying magnetic field will be re- 
placed by J2T=i Y.T=r 9{r, s). Here, g{r, s) is a generic 
function, that depends on the quantum numbers r and 
s. According to the above descriptions, the lowest en- 
ergy level is then characterized by ^^^igij", s — r), 
or equivalently by X^^i ff(''' "- = 0). 

Following the description presented at the end of 
the previous section, the fermion propagator of an 



electron in an exponentially decaying magnetic field 
is given by (see also jlj]) 



G{p,p')^{i:{p)^{p')) 



HY. 



(27r) 



Ep(p) 



{l-Pn 



m) 



rEp(p'), 



(n.39) 



where Ep(p) and p^. are given in (jll.3ip - pi.32p and 
pi.37p . respectively. Using pi.2p and the closure of 
Ep{p) from pi.33l) . it is easy to verify 

{^■n-m)G{p-p')=5^ip-p'), (11.40) 

where S-^ip^ p') = a6{t-t')6{u~u')6{y ~y')6{z~ z'). 
Here, as in the case of constant magnetic fields, since 
we are working with cBq > 0, the spin orientation in 
the lowest energy level is fixed to be positive, and the 
negative spin electrons contribute only to higher en- 
ergy levels. These facts are explicitly implemented in 
the choice 7i = s — r = 0,l,2,--- for positive spin so- 
lution F^^iu) and n = s — r = l,2,3,--- for negative 
spin solution F^Liix) in (III.29P and (III.30p . respec- 
tively.^ 



III. FERMION CONDENSATE {il^i,) IN 
EXTERNAL MAGNETIC FIELDS 

In this section, we will determine the fermion conden- 
sate (V'V') in the presence of constant and exponen- 
tially decaying magnetic fields. The fermion conden- 
sate is defined by [20| 



('0'i/') = — lim tr[Gi?(a; — x')] 



(IH.l) 



where Gf{x — x') is the fermion propagator in the 
presence of an external magnetic field, in general. In 
[5|, Schwinger's proper-time formalism is used to de- 
termine (tpip) in 8- constant magnetic field in the LLL 
approximation. In this section, we will use the Ri- 
tus method, and in particular, our results from the 
previous section to determine (ipip) first in a constant 
and then in an exponentially decaying magnetic field, 
in the LLL approximation, at finite temperature and 
chemical potential. 

A. (tptp) in constant magnetic fields 

To determine the fermion condensate in a constant 
magnetic field, let us replace G{x — x') from pi.lQp 



^ In other words, -F^^(«) is undefined and is to be neglected, 
whenever it appears in a computation. 



in piLip . Performing the traces over 7-matrices, we 
arrive at 



dpodp2dp3 



1 



x{[/;(.x)]^ + [/;_i(x)]n, (111.2) 

where according to (|TlT7|, the cUp for cBq > is de- 
fined by ujp = 2eBop+p^+m^ . In the LLL approxima- 
tion, we shall set p = 0. Using the definition of fp{x) 
from (III.16P , and after performing the integration over 

P2, 



dp2 

2tt 



Jo^i^)? 



e-Bp 

2tt '■ 



we arrive at 



imeBo f dpodp; 



2\' 



(III.3) 



(IIL4) 



where uJq 



pI 



TT J (27r)2 {j>l-uj^i 
1} . The dimensional reduction into 



two longitudinal dimensions p|| = (^07^3)1 which oc- 
curs whenever the LLL approximation is used [5| , can 
be observed in pn.4[) . To evaluate (|IIL4p . we use 



d^PlI 1 

(2^ (p2 - m2) 






(IIL5) 



which yields the fermion condensate in the presence of 
a constant magnetic field in the lowest Landau level, 



(V5^)x 



meBo 

47r2 



■In 



— ]+Oim). (111.6) 



Here, A is an appropriate momentum cutoff. This 
result coincides exactly with the result presented in 
[5[, where the method of proper-time is used. Note 
that although no integration over the coordinates x = 
(i, X, y, z) is performed here, the condensate (|IIL6p is 
constant in x. This is mainly because of the special 
form of /o (a;) from PL16I) and arises from the in- 
tegration over p2 in (jin.3|) . In Sec. IIIIBl we will 



show that in the presence of an exponentially decay- 
ing magnetic field, the condensate (ijjip) in the LLL 
depends nontrivially on x. Let us also notice that in 
massless QED, the mechanism of magnetic catalysis 
[J, |5| is made responsible for dynamically breaking of 
chiral symmetry and the generation of a finite dynam- 
ical mass in the LLL. Moreover, as it is shown in 2l|, 
apart from a dynamical mass induced by the external 
magnetic field in the LLL, an anomalous magnetic mo- 
ment is also induced, so that the dynamical mass in 
the LLL is to be replaced by the rest energy _Eo , which 
is a combination of the induced dynamical mass and 
the anomalous magnetic moment in the LLL. In this 
paper, we do not consider the effects of dynamically 
induced anomalous magnetic moment on the dynam- 
ically generated fermion mass. The latter arises ei- 
ther as a solution of an appropriate Schwinger-Dyson 
equation in a ladder approximation or by solving a 
corresponding gap equation to an appropriate effec- 
tive potential in the mean field approximation. In 
what follows, we will compute (V''0)lll from pil.6p 
at finite temperature T and finite density fi. To do 
this we use the standard imaginary time formalism, 
and replace po in pil.4p by ip^ — i{uJe, + J/i), where 
uji — 7r(2£ -I- 1)T is the Matsubara frequency labeled 
by £. Thus the integration over po is to be replaced 

by 



(ip4 1 Y^ 



(111.7) 



where P ^ = T. The fermion condensate (V'V')lll ^^ 
finite T and /i is therefore given by 



(#)^.. = - 



meBi 







7r/3 



E 



dp3 



1 



2^ (pI+pI 



(III., 



To evaluate the integration over p3 and the summation 
over £, the relation 



If. f d'^p (p2)>f 



T{l+a) ^2^^-^"+^+^^'+"^ 



(2^)^ (p2 + p2 + ^2)^ 2(47r)'^/2r(a)/3 r (f ) \P 



fc=0 



C [2{a + k - t - a) - d;^- '-^ ] + {^i ^ -fi) 



2fe 



^ I , (IIL.) 



from 23, will be used. This relation piI.Qp is orig- is the Hurwitz zeta function defined by C(s; a) ~ 
inally derived in 23| for /i := 0. In piI.Qp . C{s;a) J^kLoi^ + o)^'*, where any term with fc + a = is 



excluded. In an evaluation up to 0((m/3)''), we get 
finally 






IE 



^ 



1 z^ 

2 27r 



■i, 



1 i/3/i 

2 " "27 



^©((m/S) 



<i^-t)-^N-l^ 



(III.IO) 



where 7b ^ 0.557 is the Euler-Mascheroni number 
and il^iz) the polygamma function, defined by tpiz) = 
^lnr(2;). It arises by an appropriate regularization 
of C(l;a) using 



lim C(l + e; a) = lim "0(0) 



(III. 11) 



In pil.iop . the divergent 1/e term is neglected, and 
the identity 'tj){l/2) — —^e — 2 In 2 is used. 



B. {'4'^) in exponentially decaying magnetic 
fields 



Note that in (|III.13p . the same dimensional reduction 
to two dimensions, as in the case of constant mag- 
netic field, occurs. Evaluating the pn-integration us- 
ing pil.sp . we get 



(Vi0), 



2^ 
+0{m). 



(coshu — 1) In 



(III.14) 



In contrast to pil.6|) . the fermion condensate in an ex- 
ponentially decaying magnetic field, piI.14|) . depends 
on u = ^e"^** with 77 = x/£b, and £3 = (ei?o)~^^^, 
the magnetic length corresponding to Bq, the value 
of the magnetic field at a; = 0. Following the same 
method to introduce finite temperature and chemical 
potential as in the Sec. IIII Al we arrive at 






e""(cosh(u) - 1)<^ 21n — ^ 



-IE - 
(171/3)' 



C 3 



47r 

2 27r 

1 z^ 

2 27r 



+ tlj 



1 



2 
+ C(3; 



i/3/i^ ' 




271 J 




1 iPfi\ ' 

2 27r / 



-0((m/3)4). 



(III.15) 



The chiral condensate in an exponentially decaying 
magnetic field is given by pil.ll) . where Sf{x — x') is 
to be replaced by G{p — p') from (jTL39|. Performing 
the trace of Dirac matrices, we get first 



In Sec. |Vl we will compare the chiral conden- 
sates plLel) and piI.14|) . as well as pil.iop and 
piI.15[) . arising from constant and exponentially de- 
caying magnetic fields, respectively. 



00 00 



(V^V) = -2iemeBoY.Y. j 



dpodps 



(27r)2 (pI-u^D 

x{[F:f{u)r + [F:tAu)?}, (111.12) 



r— 1 s—r ' 



where, according to pi.37p , 



OJZ 



S,^eBoK+pl+m'^ and 



n = s^ — T^ . The expression on the r.h.s. of piI.12[) 
includes the contributions of all energy levels, denoted 
by r and s, and the terms proportional to {F'^Y ^^d 
(-^ri-i)^ correspond to the contributions of electrons 
with positive and negative spins, respectively. In the 
presence of strong magnetic fields, assuming that the 
dynamics of the system is solely determined by the 
LLL, the fermion condensate is given by 



W)lll = -2ie2^eBoe-"[cosh(u)-l] 



(27r)2 {pl-u^ir 



(III.13) 



in contrast to (IIII. 41) . The u-dependent factor behind 
the integral arises by setting n = in piI.12[) and 



using £q'' — 1 and 



Y}Fl^{u)f = e-'-Y. 



-[ (2r)! 



e " (coshu — 1) , 



IV. LOCAL ELECTRIC CURRENT 

CORRELATION FUNCTION IN EXTERNAL 

MAGNETIC FIELDS 

In this section, we will determine the local electric cur- 
rent correlation function x i^) in homogeneous and 
inhomogeneous magnetic fields. It is defined generi- 
cally by 



X«(x) 



lim 



,x. 



(i) 



(x,x') 



(IV.l) 



with 



Xn\x,x) 

= fvxVx'tT[fG^p\x,x')YGPix\x)], (IV.2) 

where i — 1,2,3 are the spatial dimensions, and the 
sub- and superscripts n on Xn and Gp indicate the 
contribution of each energy (Landau) level to the lo- 
cal electric current correlation function x (^) and 
the fermion propagator Gp, respectively. Moreover, 
Vx = dtdydz and x = [t, x, y, z). Up to an integration 



over x,^ the above definition is consistent with the def- 
inition of electric current susceptibility x presented 
in |l8|. In what follows, we will replace Gp{x — x') 
in p\r2|) by (|IlT9]) and (|lL39l) to determine the local 
electric current correlation function in constant and 
exponentially decaying magnetic fields, respectively. 



Note that in general the contribution in the transverse 
directions i = 1,2 are not equal. However, in the 
presence of strong magnetic field, where the dynamics 
of the system is reduced to the lowest Landau level, 
p = 0, and we have therefore 



AA„«=AAf =0, 



(IV.. 



A. x'''(^) i*! constant magnetic fields 

For a constant magnetic field, the contribution of 
Xp \x, x') corresponding to each Landau level, labeled 
by p, is given by 



X^{x,x') = / VxVx'tY[-i'Gp{x,x')^'Gp{x',x)\, 



(IV.3) 



where Vx — dtdydz. Using pi. 191) and plugging 



Gp{x,x') = Vp e'P-(*-*')Pp(a;) 



J -Pp-m 



-Pp{^'), 



(IV.4) 



with Pp{x) from (JIlTTI) and Pp from (IIlTtI) in (JTOI) . 
we arrive after integrating over coordinate and mo- 
mentum variables, and taking the limit x' — >■ x, first 
at 



xi'H^)^ \imxi'H^,x') 
Vp 



,LyLz 

'' T 



{Pi 



,2^2 



tilYPpix) {'ypp + m) Pp{x) 



X 7*Pp(x) (7 • Pp + m) Pp{x)] 



(IV.5) 



Here, Ly and Lz are constant lengths of our probe in 
second and third spatial directions. They arise from 
the integration over y and z, respectively. The tem- 
perature T in pV.sp arises from the integration over 
the compactified imaginary time t = it (z [0, /3], where 
/3 = T~^. Performing then the trace over Dirac ma- 
trices, and summing over all Landau levels, we get 



x'^'Hx) =2i 



00 



L,,L 



T 



dpQdp2dp3 



p^,. i2nr [Pl-^l? 



^^ , (IV.6) 



where uj'i = p\+Pz + rn? with p\ — 2eBf)p for eB^ > 
[see (|lL17l) ]. and 

AAp(i) = 2ipi-u:i)[f+ix)nf+_,ix)r, 

AAf = 2ipl-u;l+pl)[f+ix)f[f+_,{x)f, 

U;,'^ = [pl-cal + 2{pl~plm+{xr 

+ [pI -C.I + 2ipl +plm+_,{x)]\ (IV.7) 



The X direction is specified, because of the specific Landau 
gauge, which fixes A^ as is described in Sec. Ill Al 



which lead to vanishing electric correlation functions 
in the transverse directions, i.e. Xo = 0- ^s concerns 
the electric current correlation function in direction 
parallel to the external magnetic field, Xq, plugging 

Ml,'^=ipl-u:l + 2pl)[f+ix)]\ (IV.9) 



from pV.7|) in (IIV.6p . using the definition of f^ix) 
from pi.l6p . and integrating over p2, 



/tW(.r^(f 



3/2 



(IV.IO) 



we arrive at 



Xo = 



2lLyLz 

7^~ 



eSo 

27r 



dpodps, ipl - wg + 2pl) 



(2^)2 (pl~ul) 



2\2 



(IV.ll) 



Introducing the temperature as in the previous sec- 
tion, and using (IIII.9I) to perform the sum over I, la- 
beling the Matsubara frequencies in 



1 v^ [ dp3. 


1 

Xpl+pl + m-^) 


•^pI 


{pl+pl+m^Y 


1 




(IV.] 



we arrive finally at the components of the electric cor- 
relation function in the transverse and longitudinal 
directions with respect to a constant and strong mag- 
netic field in the LLL approximation 



X^ 



Xo 



0, 

LyLz 

2tiT 



2ti 



3/2 



(IV.13) 



Note that whereas Xo" — 0, the non-vanishing Xo does 
not depend on x, although no integration over x is per- 
formed in the definition pV.3p of Xp {x) from (IIV.SP . 
To arrive at pV.13p . the resuh from pV.12p is used. 
In the chiral limit ?77, -^ 0, both integrals appearing on 
the l.h.s. of pV.12p are infrared (IR) divergent, and 
the (bare) fermion mass to, plays the role of an IR 
regulator. Using the relation pil.9p to expand these 
integrals in a high-temperature expansion, the IR di- 
vergences of the integrals appearing in pV.12p cancel 



and we arc left with an exact, ni independent solu- 
tion for Xg, which does not also depend on the chem- 
ical potential /i. Note that the above result (|IV.13p 
suggests that for cBq <C T"^ 



and LyLz 



< 



£% with 



£b — (eSo)^^^^, apart from x^, Xo also vanishes. 

The electric susceptibility of a massless magne- 
tized QED at zero temperature is recently studied in 
24| . According to the mechanism of magnetic cataly- 
sis [J, la], in massless QED, the exact chiral symmetry 
of the original QED Lagrangian is dynamically bro- 
ken by the external magnetic field. In [24| , it is shown 
that the electric susceptibility in the chirally broken 
phase at zero temperature is independent of the ap- 
plied magnetic field. It would be interesting to extend 
the results presented in [2J| to finite temperature and 
study the temperature dependence of electric suscep- 
tibility in chirally broken and symmetric phases of a 
magnetized QED. This is indeed an open question and 
we hope to come back to this point in the future. 



B. X (w) in exponentially decaying magnetic 
fields 

In the case of exponentially decaying magnetic fields, 
the dimensionless coordinate u plays the same rule as 
X in the previous section. Thus, defining 

xi') {u, u') = J VxVx'ty[YG.{p, p'h'G^p' , p)] , 

(IV. 14) 
in analogy to (jIV.Sp , using (JII.39P and plugging 



For different directions, i — 1,2,3, the nominator Mr,l 
is given by 

A/;W = 2{pl-ul)[Fl^{u)nFlL,{u)]\ 

Af^ = [pl-ujl + 2{pl-pl)][F?:{ur 

+ [pI -0^1 + 2{pI+pI)][F^:_,{u)]\ (IV.18) 

where F^^{u) is defined in pi.29p . In the LLL, where, 
according to our explanations in Sec. IIIBl s — r and 
therefore k = 0, the only contribution to x''*'' arises 
from the positive spin particles. We get therefore 



AC« =A/-il. = 0, 



(IV.19) 



leading to vanishing transverse components of the 
electric current correlation function in the LLL ap- 



proximation, i.e. Xl 



0. The same effect is also 



observed in the case of constant magnetic fields in see 
Sec. IIV Al Moreover we get, 

<1. = {pI u^l + 2pl)[F^^{u)]\ (IV.20) 
Plugging PV.20I) in (IIV.17I) . we arrive first at 



xL.(«) 



2ia^LyL 



y-^^- \^lTp2r 



T 



r=l 



F^^{u)r 



dpodpz {pI -ujI + 2pl) 



(27r)2 (p^-u;(j) 



2^,2 



(IV.21) 



where || denotes the third (longitudinal) component of 
X^*-* in the LLL approximation. In (IIV.21[) . as in Sec. 
IIIIB| the summation over r can be performed using 



y^Prfin) 



P^^iu'), (IV.15) 



J ■ Pk- m 

with P^''{x) from pL32]) and p^ from (|IL37)) in 
(jIV.14p . we arrive after taking the limit m' — ?> m at 



Xi^n«)= limxl^H^,"') 

u'—^u 



dpodp-i 



1 



ia LyLz 
^^T^J ~{2^Jpl-u:ir 
xtr[YP^^iu){^-p, + m)P^^iu') 
•KfP^^iu') (7 . p, -f m) Pl^iu)], (IV. 16) 

where a = £,^/eBo, n = s—r, and, according to (III.37p . 
^K = P2+P3 + "^^ with P2 = i^eBfjK and k = s^ — r^. 
After performing the trace over Dirac matrices, the 
local electric current correlation function including the 
contributions of all energy levels reads 



00 00 



x«(")-^^^EE 



T 



r— 1 s^r 



dpodps 



j4!I 



{2ny (p2-^2) 



2^2 ■ 



(IV.17) 



r=l 



[Ftiu)Y 



00 / o^ 



^\{2r)\ 

r— 1 ^ ^ ' 

[(/o(2u) + Jo(2u))-2],(IV.22) 



and C^ = 1. Here, l{){z) and Jq{z) are zeroth order 
modified Bessel functions / and J . Plugging (IIV.22P 
in pV.2ip . we arrive at 



Xlll(w) 



ia LvL 



T 



y^z^-2u 



[{Io{2u) + J„{2u)) ~ 2] 



dpodps (pI -uj^ + 2pI) 
' (2^)2 {pl-.^,r ■ ^^ '^ 



Introducing the temperature and finite density, using 
the method introduced in the previous section, and us- 
ing (|IV.12p to perform the integration over p^ and the 
sum over the Matsubara frequencies, the transverse 
and longitudinal components of the local electric cur- 
rent correlation function in the LLL approximation 
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read 



Xlll 



xL(w) 



0, 



4:TTT 



xe-2" [/o(2m) + M2u) - 2] . (IV.24) 



In the next section, we will compare local electric cur- 
rent correlation functions PV.ISP and (|IV.24p arising 
from constant and exponentially decaying magnetic 
fields, and discuss the remarkable property of (jIV.24[) 
in the limit ^ — > 0. 



V. DISCUSSIONS 

In the previous sections, the chiral condensate (V'V') 
and local electric current correlation function x ^I'e 
computed in the presence of constant and exponen- 
tially decaying magnetic fields in the LLL approxima- 
tion. As it turns out the transverse (z = 1,2) compo- 
nents of X vanish in this approximation. Moreover, 
whereas in the presence of constant and strong mag- 
netic fields (V'V')lll from (|nL6l) and xl! from (|IV.13P 
are constant (as a function of coordinates), they de- 
pend, in the presence of exponentially decaying mag- 
netic fields, on a variable u defined hy u — -pje~^^ . 
Here, rj is proportional to the coordinate x, and is 
given explicitly by 77 = x/Hb, where the magnetic 
length Hb = (ei?o)^^^^- To give an explicit exam- 
ple on 77, let us take cBq — 15to^ ^ 0.3 GeV^, with 
the pion mass TTi^r ^ 140 MeV. This corresponds to 
a magnetic field Bq ^ b x 10^^ GauB,^ which is the 
typical magnetic field produced in the early stages of 
heavy ion collisions at LHC 26[ or exists in the in- 



terior of compact stars [27[. In this case, the corre- 
sponding magnetic length is given by is ~ 0.4 fm.® 
Taking rj = A, for instance, would mean a distance 
X = Mb ~ 1.6 fm from the origin at x = 0, etc. Note 
that for the above mentioned LLL approximation to 
be reliable, we have to consider only small damping 
parameter ^ and consider the probe in small rj from 
the origin a; = with B{x — Q) — Bq. 

In this section, we will compare the values of 
(^V)lll from (|nL6|) with (V5V)lll from piI.14p . as 
well as Xo^ from (|IV.13|) with xlli. from (|IV.24|) . in 
order to explore the effect of inhomogeneity of the 
external magnetic field, especially once quantum cor- 
rections are taken into account. To do this, let us first 



''' We are working in the units where eB ■■ 
to Bo ~ 1-7 X lO^o Gaufi [H. 



1 GeV corresponds 



As it turns out, the magnetic length Cb for eBo = 1 GeV is 
given by Ig ~ 0.63 fm. 




FIG. 1: The ratio ^(77) = (i/'V'>lll/(i/'V'>£ll is plotted as 
a function of r; = x/(.b with £s = (eBo)"^ for different 
damping parameters ^ = 0.5, 1, 1.26, 2.5. The condensates 
(7/'7/))lll and (i/iV')£LL are defined in (1111.14^ and pil.6|l . 
respectively. For the damping factor ^ ^ 1.26, the ratio 
C(_(ri) is maximized at r; = (or equivalently a; = 0). 



define the ratio 



Q(^) 



_ (7^^), 



{i,^)[ 



= 27r^2g-n (cosh(M) 



1), (v.i) 



with u = 4ie^^^ and rj = xJIb- 

In Fig. [U C^ijj) is plotted for different damping pa- 
rameters ^ = 0.5, 1, 1.26, 2.5 as a function of 77 = x/Ib- 
Depending on the damping parameter S,, the conden- 
sate ('!/'7/')lll arising from an exponentially decaying 
magnetic field is up to a factor 2.5 greater than the 
condensate (V''0)lll m the presence of a uniform mag- 
netic field. For ^ ~ 1.26 the maximum value of the 
condensate (V'V')lll arises at 77 = or equivalently at 
a; = 0. Let us notice that this interesting observation 
is indeed in contradiction to our prior expectation, 
according to which we expect Q(7y) = 1 at 77 = (or 
equivalently at x — 0), because B{x — 0) — Bq is 
constant. This indicates the singular nature oi rj — 
(or equivalently a; = 0), once quantum fiuctuations 
produce a nonvanishing chiral condensate (tpip)- 

The singular nature of Q (77) for ^ = is explored 
in Fig. [21 where the same ratio ^(77) is plotted as 
a function of the damping factor ^ for 7/ = 0, 0.5, 1, 2 
(or equivalently for x = 0,£b/'2,^b,'2^b)- For 7/ ~ 0, 
Q(0) is maximized at ^ ^ 1.26. For 7; 7^ 0, how- 
ever, the maxima of ^(77) are shifted to smaller val- 
ues of £,. At ^ = 0, the condensate {^p'lp)^^^, arising 
from an exponentially decaying magnetic field, van- 
ishes for all demonstrated values of 77. This is again 
in contradiction to our prior expectation, according to 
which we expect ^(77) = 1 for ^ = 0. This is because 
for '^ = the exponentially decaying magnetic field 



11 



2.5 


/N 


_i 


/ \ 


"i^.O 


/ \ 


^ 


/ \ 


'?1.5 


K \r] = 


Jl 


9 \ 


_l 


uv » \ 


^1.0 


A \ \ 


1^ 


ft \ \ Nv 


^0.5 


/ \ \ \ ^^^^ 




/ \ \ ^ ^^*"^^*-—^ 




/ \ \ ■* ''^—-^~ _ 


A n 





10 



^ 



FIG. 2: The ratio ^(ry) = (iAV>lll/{^iA)lll from (|Vl|) 
is plotted as a function of tiie damping factor ^ for 77 = 
0,0.5,1,2 (or equivalently for x = 0,£B/2,£B,2eB witii 
fs = (e_Bo)~^ ) from above to below. For r; — 0, the 
^(77) is maximized at ^ ~ 1.26. For i] j^ (or equivalently 
a; 7^ 0), the maxima are shifted to smaller values of ^. For 
all values of r?, the condensate {'tptp) 1,1^1,, arising from an 
exponentially decaying magnetic field, vanishes at ^ = 0. 




FIG. 3: The ratio B{x)/Bo = e"''« with 77 = x/£b and 
£b = {eBo)~^'^ , is plotted as a function of the damping 
factor rj for ^ = 0, 0.5, 1, 2 from above to below. 



B{x) = Bqc^^'' becomes constant (see Fig. [3]). 

The same singular behavior at 77 = and ^ = is 
also observed when we repeat the above analysis for 
the ratio 



xdv) 



xL(«) 



Xo 
= V2^e^e-'"[/o(2u) + Jo{2u) ~ 2], (V.2) 



where xU from (|IV.13p and Xi.ll(") from (|IV.24I) are 
local electric current correlation functions in the pres- 
ence of a uniform and an exponentially decaying mag- 
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FIG. 4: The ratio xdv) = Xlll(")/Xo from (EU is plot- 
ted as a function of 77 = x/£b with £b ~ (e_Bo)~^ for 
different damping parameters ^ = 0.1,0.2, ■■• ,1,1.26,2. 
For ^ '^ 1.26, the ratio X(.{v) is maximized at 77 = 0. 



^e" 



^ri 



netic fields, respectively. Here, as before, u 
and r] = x/£b, with £3 = (eBo)"^/^ 

In Fig. m the ratio x^(?7) is plotted for ^ — 
0.3,0.4, ••• , 2 as a function of 77. The fact that 
X^iv) 7^ 1 for ?/ = is again in contradiction to 
our prior expectation: We know that the x-dependent 
magnetic field B{x) = B^e^^^ becomes constant for 
77 = (or equivalently a; = 0), and therefore the prior 
expectation is that Xi iv) = 1 as well as Q (77) = 1 for 
?7 = 0. The above observation demonstrates again the 
singular nature of the limit x — > 0, once the quantum 
effects are to be considered. Another remarkable point 
here is that, according to Fig. [51 for 77 = (or equiv- 
alently X — 0) and in the interval 0.5 < ^ < 4, the 
ratio C^{ri) > 1, whereas according to Fig. |3](see also 
Fig. [SJ, Xii'u) ^ 1 for 77 = and in the same interval 
0-5 ^ C ^ 4. In other words, for r; = and in the 
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interval 0.5 < (, < 4,^ the value of the electric current 
correlation function Xo arising from a constant mag- 
netic field is always larger than XlllC^) arising from 
an exponentially decaying magnetic field, whereas in 
the same interval of ^, (V'V')lll arising from a constant 
magnetic field is always smaller than ('!/''0)lll arising 
from an exponentially decaying magnetic field. 

In Fig. [5l the same ratio x^iv) from (jV.2p is plot- 
ted as a function of the damping factor ^ for fixed 
7; = 0, 0.5, 1, 2 (or equivalcntly for x = 0, £3/2, Ib, 2^s 
with (.B = (eSo)^^^^) from above to below. As it 
turns out the value of Xlll(u) for ^ — > vanishes, in 
contrast to our expectation. 




FIG. 5: The ratio Xii'n) from (|V.2p . is plotted as a func- 
tion of the damping factor ^ for r) = 0, 0.5, 1,2 (or equiva- 
lently for x = 0,^s/2,^s, 2^s with £b = (eBo)-^''^) from 
above to below. For r; = 0, Xiiv) is maximized at ^ ~ 1.26. 
For rj ^ (or equivalcntly a; 7^ 0), the maxima are shifted 
to smaller values of ^. For all values of rj, Xlll (w) vanishes 
at ^ = 0. 
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FIG. 6: The values of ^ that maximize the ratios C^{ri) 
from (|V.H) [blue circles] and Xiiv) from (|V.2|) [green rect- 
angles] are plotted as a function of 77. They almost overlap. 



ters ^ = 0.1,--- ,10. Each dot in Figs. [7^-[71; indicates 
a value of ^ = 0.1,- •• ,10 with A^ = 0.05. As it 
turns out, Xiiv) increases with increasing Q(?7) for a 
certain numbers of ^ (see Fig. [7^). The maximum 
value of ^ for which this is true, is different for dif- 
ferent 77. Let us denote this specific £, with f^,. For 
ri = 0,0.2,0.5,1, C*'s are ^^ - 1.3,1.1,0.95,0.3, re- 
spectively. For ^ > ^i,, however, x^iv) decreases with 
decreasing C^{r]) (see Fig. [7)d). In Fig. ^, the two 
curves in Fig. [7^ and [7b are reassembled and demon- 
strate a hysteresis-like curve with increasing ^. How- 
ever, since the LLL approximation is only valid for 
small ^ and rj, we shall limit us to the values of ^ < ^^ 
for which Xiiv) increases with increasing Q(r?) (Fig. 



m- 



The second observation from Fig. [5] is that the 
maxima of Xiiv) are shifted to smaller values of ^ 
for increasing 77. In Fig. IH the values of £, that maxi- 
mize the ratios defined in (jV.ll) [blue circles] and (|V.2[) 
[green rectangles] are plotted as a function of 77. They 
are almost the same. This indicates a certain rela- 
tion between the formation of chiral condensates and 
the value of electric current correlation function in a 
system including relativistic fermions. 

This relation is demonstrated in Figs. [7^-[71;, where 
X5(77) from (|V.2p is plotted as a function of ^(77) 
from (|V.1[) for fixed 77 = 0, 0.2, 0.5, 1 (or equivalently, 
X — 0,lB/5,iB/2,iB) and various damping parame- 



^ As before mentioned, the LLL approximation is only reliable 
for small values of ^ and rj. 



VI. CONCLUDING REMARKS 

In the present paper, the effect of exponentially decay- 
ing magnetic fields on the dynamics of Dirac fermions 
is explored in detail. Using the Ritus eigenfunction 
method, we have first determined the energy spec- 
trum of fermions in this non-uniform magnetic field 
and compared it with the energy spectrum of rela- 
tivistic fermions in a constant magnetic field. We 
have then computed the chiral condensate {^p^p) and 
local electric current correlation function x in the 
presence of strong uniform and non-uniform magnetic 
fields in the LLL approximation. In non-uniform mag- 
netic fields, {ipip) and x^*'' depend on a dimensionless 
variable u = ^e~'>^'^^, which is a nontrivial function 
of the magnetic field's damping parameter ^, the co- 
ordinate X and the magnetic length (.b = {eBo)~^/^. 
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FIG. 7: The ratio Xe('?) from (|V.2|) is plotted as a function of the ratio C^{ri) from HV.ll) for fixed tj = 0,0.2,0.5, 1 
(or equivalently, x = Q,(.b/^,^b/'2,(.b) and various damping parameters ^ = 0.1, ■ • • ,10. Each dot in the above curves 
indicates a value of ^ = 0.1, • ■ • ,10 with A^ = 0.05. For ^ smaller than a certain ^j,, xsl'?) increases with increasing 
C(_{r]) (panel a). The value of ^* is different for different 77 (see the main text). For ^ > ^*, however, Xiiv) decreases with 
decreasing C^{r]) (panel b). The two curves in panels a and b are reassembled in panel c. In the LLL approximation, for 
small values of ^ and rj, only the behavior demonstrated in panel a is relevant. 



In the LLL approximation, the transverse components 



of the electric susceptibility, Yllli ^ = 1; 2, vanish and 
its longitudinal component Xlll depends also on the 
variable u, only in the presence of exponentially de- 
caying magnetic fields. We have shown that the lim- 
its ^ —> as well as x — > are singular. In these 
limits, the x-dependent magnetic field becomes con- 
stant, and therefore ('0'0)lll as well as Xlll arising 
from the a;-dependent magnetic field are expected to 
have the same value as in a constant magnetic field. 
But, as it turns out this is not the case. This remark- 
able behavior of {ipip) 1.1.1. as well as Xlll at ^ — >■ 
and/or a; — ?> is discussed in detail in Sec. |Vl Let us 
notice that, mathematically, there is a difference be- 
tween taking the limit of a quantity to zero and setting 
it exactly equal to zero. This difference is indeed re- 
sponsible for the singular behavior of the limits 1^ — > 
and/or x — >■ 0. When we are looking for the solution of 



the differential equation (jll.221) for non-uniform mag- 
netic fields, for instance, the parameter ^ appearing 
in u = ^e^''^/^'^ can be very small (^ — >■ 0), but 
it cannot be set exactly equal to zero. When £, — 0, 
the magnetic field becomes constant, and the differ- 
ential equation (jll.22|) . leading to the solution (11.31), 
including associated Laguerre polynomials, has to be 
replaced by the differential equation (JII.ISP for con- 
stant magnetic fields, leading to Hermite polynomi- 
als. There is no way to reproduce the Hermite solu- 
tion of the latter case from the associated Laguerre 
solution of the former case, by taking ^ — >■ 0. This is 
why that although at the classical level the uniform 
magnetic field is reproduced in these limits from the 
non-uniform magnetic field, neither the energy spec- 
trum nor the quantum corrections to (?/"/') lll as well 
as Xlll of the constant magnetic fields can be repro- 
duced from their values in the non-uniform magnetic 
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n ~ n > r 



n > n and r > r 
n — n = r ~ r n~n<r— r<n n < r 



n > n' and r' < r 
r — r < n r — r > n 



h 



h 



h 



h 



TABLE I: The integrals Jj, z = 1, ■ • • , 6 are presented in (|A.3|) and (|A.4|) . 



field by taking the limit ^ — >■ 0. 

We have also shown that for small values of x/^b 
and ^, where we believe that the LLL approximation 
is reliable, the ratio xcl'?) from (jV.2p increases with 
increasing ratio C^{r]) from (jVTp up to a maximum 
value of ^, denoted by ^*. As it turns out, the value 
of ^^, depends on -q (see Sec. IVT) . 

Let us notice that electric susceptibility (electric 
current correlation function) of a three-flavor color su- 
perconductor is recently computed in a strong homo- 
geneous magnetic field in 28] , and the magnetoelectric 
effect in a strongly magnetized quark matter is stud- 
ied in detail. The present work is a second nontriv- 
ial example on the effect of external magnetic fields, 
in general, and spatially decaying magnetic fields, in 
particular, on electric current correlation function of a 
system containing relativistic Dirac fermions. Apart 
from applications in condensed matter physics, the 
results of this paper may also be relevant in the con- 
text of heavy ion collisions. Here, recent experimental 
activities at RHIC indicate the production of intense 
magnetic fields in the early stage of non-central heavy 
ion collisions 26|. Depending on the initial condi- 
tions, e.g. the energy of colliding nucleons and the 
corresponding impact parameters, the spatially vary- 
ing magnetic fields are estimated to be in the order of 
B ^ 10^^ — 10^^ Gaufi, decaying very fast, within few 
femtometers to S ~ 10^^ — 10^'' Gauf5, and vanishing 



J 



during the hadronization process [26 



VII. ACKNOWLEDGMENTS 

The authors acknowledge valuable discussions with 
Sh. Fayazbakhsh, S. A. Jafari and A. Vaezi. 



Appendix A: Orthonormality relations of 
associated Laguerre polynomials 

In this appendix, we present a number of useful or- 
thonormality relations of associated Laguerre polyno- 
mials. To do this, let us define the integral 



/: 



2r' / 



du e-"u''+'' ct;{u)ci': (u). (a.i) 



For r = r' , the standard orthonormality relation of 
associated Laguerre polynomials is given by 



/: 



(n-t-2r)! 



(A.2) 



But, in the present paper, we have to consider the 
cases, where in general r ^ r' . To determine /^'^/, we 
have to distinguish several cases, which are summa- 
rized in Tables HI and HIl 



n < n and r' < r 
n — n > r — r n — n = r — r n— n<r~r<.n n < r ~ r 



n < n' and r' > r 
r — r < n r — r > n 



h 



is 



J9 



'10 



111 



/l2 



TABLE II: The integrals Ii,i^ 7,- ■ ■ ,12 are presented in (|A.5|) and (|A.6|) . 



The integrals /j , i = 1 , • • • , 4 for n > n' and r' > r are given by 

h = 0, 

(-l)'^'-'' {r' +r + n'y. 



h = 



_ ^'■' '■y^" "'^ (-1)"-"'+^ (/ - r){r' +r + l- n')\ {r' - r + I - 1)! 
■■^ ~ £! (n' - €)! {n-n' + (.)\ {r' -r-n + n' - £)! ' 



1=0 



^^ = E 



«=0 



(-l)"-"'+^(r' - r){r' +r + £- n')l{r' - r + £ - 1)1 
£\ {n' - l)\ (n - n' + i)\ (r' - r - n + n' - t}\ 



(A.3) 
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For n > n' and r' < r, we have 



E 



(-l)^(r - r'){r' + r + n' ~ £)\{r - r' + n ~ n' + £ ~ 1)1 
i\ {n' - e)l{n -n' + ^)!(r - r' - 1)1 ' 



^ (-l)^(r - r'){r' + r + n' - £)!(r - r' + n - n' + £ - 1)1 

h = 2^ 

t=0 



£\ (ri! - £)l{n - n' + £)\{r - r' - £)] 
The integrals Ii,i = 7,- ■ ■ ,10 for n < n' and r' < r are given by 



h = 0, 



(-!)'■-'' (r' + r + n)! 



''■ '''y^"' ") (_!)»'+£-» (^ _ /)(^ + / - n + £)! (r - r' + £ - 1)! 
"^ ~ £! (n - £)! {n' -n + £)l {r - r' - n' + n - £)! ' 



e=o 



e=o 
For 71 < n' and r' > r, we get 



(-l)"'-"+^(r - r')(r + / + £ - n)!(r - r' + £ - 1)! 
£\ (n - i)\ (n' -n + £)! {r - r' - n' + n - £)[ 



"^ {~l)^{r' - r){r' + r + n- i)\{r' - r + n' - n + £ - 1)1 

hi - 2^ 

1=0 



£1 {n - e)\{n' -n + e)\{r' - r - £)l 



hi = 



E 



(-l)^(r' - r){r' + r + n- <')!(r' - r + n' - n + £ - 1)! 
£\ {n - £)\{n' -n + £)\{r' - r - £)\ ' 



r 



(A.4) 



(A.5) 
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